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Abstract 



A relativistic extension of our pseudo-shifted ^-expansion technique is pre- 
sented to solve for the eigenvalues of Dirac and Klein-Gordon equations. Once 
more we show the numerical usefulness of its results via comparison with avail- 
able numerical integration data. 



I. INTRODUCTION 



Dirac and Klein-Gordon (KG) equation are not exactly soluble for most of Lorentz vector 
( coupled as the 0-component of the 4- vector) and/or Lorentz scalar ( added to the mass 
term) potentials[l-14]. One, therefore, has to resort to some approximation schemes [1- 
6]. Yet, in between non-numerical and purely numerical non-relativistic (Schrodinger) and 
relativistic (KG and Dirac) wave equations there exists a broad gray zone of potentials 
tractable via various systematic semi-numerical (or semi-analytical) power-series expansions. 

In numerous methodical predecessors of a subset of papers, Mustafa and co-workers [7] 
have sought a possibility in the power-law asymptotic expansions using some small param- 
eter to solve for Schrodinger equation. It has been noted that the presence of the central 
spike Id {^d + 1) / (where la = ^ + (c? — 3)/2 and the dimensions d > 2) va. the radial 
Schrodinger equation, just copies the effect of the centrifugal and/or centripetal force and 
immediately inspires the use of small shifted inverse angular momentum quantum number. 
Their PSLET ( pseudo-perturbative shifted-i? expansion technique) has provided persuasive 
numerical verifications by immediate comparison of its results with available brute force nu- 
merical data [7]. PSLET simply consists of using 1/1 as a perturbative expansion parameter, 
where T = £ — /3o, £ is a quantum number, and (3o is a suitable shift introduced to avoid the 
trivial case i = [7].pagebreak 

In this paper, we extend PSLET recipe to solve for Dirac and KG equations with Lorentz 
scalar and/or Lorentz vector radially symmetric potentials (in section 2). In section 3 we 
apply this relativistic recipe to some exactly solvable potentials ( e.g., V{r) — S{r) = —A/r, 
V{r) = —Ai/r and S{r) = 0, V{r) = and S{r) = —A2/r, and the Dirac oscillator) and 
non-exactly solvable (by PSLET) potentials ( e.g., the pure scalar linear, the funnel- shapped, 
and the power-law potentials) to study the usefulness of its numerical results. We conclude 
in section 4. 



The Dirac equation with the Lorentz scalar ( added to the mass term) and Lorentz vector 
( coupled as the 0-component of the 4- vector potential) potentials reads {mh = c=l units) 



II. PSLET RECIPE FOR DIRAC AND KG EQUATIONS 



{a.p + P[m + S{r)]} *(r) ^ {E - V{r)} *(r). 



(1) 



Which decouples into 



Ur)G{r) + 



dF{r) 
dr 



F{r) = 0, 



(2) 



2 



6 r F r - -G r 

dr r 



0. 



(3) 



where k = - + 1) for j = ^ + 1/2, k = ^ for j = £ - 1/2, and 

ei(r) = E-1/(r) - [m + S{r)], 



^2{j') = E + m — y{r); y{r) = V{r) — S{r). 

E is the relativistic energy, and G{r) and F{r) are the large and small radial components 
of the Dirac spinor, respectively. In terms of the large component G(r), equation (2) reads 



d k(k ~t~ l) 1 , I , ^ r d K-, ^ J. / \ J. / \ 



dr 



6(r) 



dr 



G{r) = 0, 



where the prime denotes d/dr. It can be shown that with the ansatz 

G{r) = $(r) exp{-p{r) / 2) ; p'(r) = y'(r) / ^2(0, 
equation (4) reads 



dr"^ 



+ ' , ' + U{r) - 6(r) $(r) = 0, 



(4) 



(5) 



(6) 



where 



2 6(r) r e2(r)^4 U2(r). 



(7) 



Obviously, equation (6) reduces to KG-equation with k + 1) = + 1), for any k, if 
Uir) is set zero. It is therefore convenient to introduce a parameter A = 0, 1 in f/(r) so 
that A = and A = 1 correspond to KG and Dirac equations respectively. Also, we shall 
be interested in the problems where the rest energy vtk? is large compared to the binding 
energy — E — mc? . This would manifest the approximation 
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1 



i2{r) Eund.^'^m-yir) 2m 



0{l/m% 



(8) 



which in turn implies 



U{r) 



A 
4m 



y [r) + — j — 

r 4m 



(9) 



For Coulombic-like potentials (i.e., a Lorentz-vector V{r) = —Ai/r and a Lorentz-scalar 
S{r) — —A2/r potentials) one may re-scale the potentials and use the substitutions 



Vr{r) = V{rf - §, 



(10) 
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(11) 



to recast equation (6) as 

, [P + r(2/3„ + l)+/34/5, + l)] 



+ r{r)+2EV{r) 



where 



r(r) = + Sr{r) + 2mS{r) + m^ + U{r), 



:i2) 



(13) 



(14) 



and is a suitable shift to be determined below. Next, we shift the origin of the coordinate 
system x = l^^'^{r — ro)/ro, where Tq is currently an arbitrary point to be determined through 
the minimization of the leading energy term below. It is therefore convenient to expand 
about X — (i.e., r = Tq) and use the following expansions 



,-2^Y,0^^nJ-n/2. ^ {-if (n+l). 



n=0 "^ o 



(15) 



72 oo rl''^V(r \ 



(16) 



Vixir)) = V c x"r"/'- c = ^''^i^o) !l 



(17) 



-1 oo 

n=-l 



(18) 



where Q is set equal to P at the end of the calculations. With the above expressions into (12), 
one may collect all x-dependent terms of order / to imply the leading-order approximation 
for the energies 



Which upon minimization, i.e., dE^'^^dro — and d^E^~^^ /drl > 0, yields 

2Q^h {To) + \lh{Tof-g{To) 

where 

h{r,) - rl \2 V{to) V {to) + r'(r„) + roV\r,f 



(19) 



(20) 



(21) 
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(22) 



and primes denote derivatives with respect to Vg. This imphes that xl ^-coefficients vanish, 
i.e.; 



Qai + rlhi + 2rlE^-^^ Ci = 0. 
Equation (12) therefore reduces to 



'2 oo 



+ ^ r„ x^r + (2 /3„ + 1) ^ a„ x'^ I- 



■n/2 



n=2 



n=0 



n=0 

9»,2 oo n+1 

^ ^ £;(n-p) 1^^^^ ^2p J-n ^ ^^^^^ ^2p+l^-(n+l/2)^ 



Q n=0 p=0 

— ^ ^ ^(n-p) ^(p) J--(n+l) 
Q n=-lp=-l 



(23) 



(24) 



where 



(25) 



One may now compare equation (24) with the Schrodinger equation for the one-dimensional 
anharmonic oscillator 



d^ 1 



+ -u;V + £o + 5(y) 



d'lp' 4 

where is constant, i^(?/) is a perturbation like term and 



n=l 



with A; = 0, 1, 2, • • • and 



/ 2r4 4 

a; = Ji2 + -f r"(r„) + — ^ £;(-!) y"(r„). 



g ^ " Q 

In a straightforward manner, one can show that 



Q 



2rl [£;(-!) - Co] 



[(2/3„ + l) + (fc + l/2)a;], 



(26) 



(27) 



(28) 



(29) 



and choose Po so that E^^^ = to obtain 



equation (24) then becomes 



n=0 

reak where 



v^'^\x)^T2x'' + {2po + l) ao, 



(30) 



$,,,(x) = 0, (31) 



(32) 



v^^\x)^T3X^ + {2Po + l) aix, 



(33) 



V 



(x) = 7;+2 + (2 /3„ + 1) a„ x^ + (3o (/3„ + 1) a„_2 x'^"'; n > 2, (34) 



/o»,2\ n+1 



(35) 



/o»,2\ n+1 



(36) 



An) ^ Ijl 

Q 



n+1 

/ p=-i 



Now we may closely follow PSLET recipe for the fc-nodal wavefunction and define 

$jk,^(x) = Fk,e{x) exp{Uk,i{x)) 

where 



oo k—1 



(37) 



(38) 



(39) 



n=0 p=0 



(40) 



n=0 



with 



n+1 

Ul:kx) = E D,,r^,kx'^-'■, Do,n,k = 0, (41) 
p=0 



n+1 

Gt^eix)-T.Cp,n,kx'^. (42) 

p=0 

Equation (31) then reads 

oo 
n=0 

-F,,,(x) [C/;,,(x) + Ul,{x)Ul,(x)] -2FlAx) Ul,{x) - fI,{x) = (43) 

where primes denote derivatives with respect to x. One may also ehminate /-dependance 
from equation (43) to obtain four exactly solvable recursive relations ( see Appendix for 
details). Once To is determined, through equation (20), one may then calculate the energy 
eigenvalues and eigenfunctions from the knowledge of Cp^n,k, Dp^n,k and A^^l in a hierarchical 
manner. 



III. ILLUSTRATIVE EXAMPLES 

In this section wc illustrate the applicability of the above relativistic PSLET recipe 
through some examples covering Dirac and KG equations. 



A. An equally mixed Coulomb potentials 

For an equally mixed Coulombic potentials, i.e. V{r) = S{r) = —A/r, U{r) in (9) 
vanishes. Consequently r(r) = —2mA/r + m^, Q = —fi? + 2mAr, uj = 2, f3o = — {k + 1), 



{£' + k + 1) + A"^ I (2m74), and the leading-order approximation reads 

2^2 



m 



1 - 



(A; + £ + 1)^ + ^2 



which is the well known exact result for the generahzed Dirac- and KG-Coulomb problems, 
where higher-order corrections vanish identically. 



B. Vector Coulomb or scalar Coulomb potential 

For V{j-) = —Ai/r and S{r) = or V{r) — and S{r) — —A2/r in KG equation one 
would obtain the well known exact results 



-1/2 
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or 



2 



1/2 



; n2^k+^- + ^{£ + l/2f + Al 



respectively. Again higher-order corrections vanish identically. 



C. Dirac oscillator 



Following the work of Romero et al [13], the Dirac oscillator [14] eigenvalue problem 
see equation (30) in [13]) reduces to 



<f ^ A(A + e/?) 



+ rrfB^r^ + + mB (e [2j + 1] - /3) 



where A = j + 1/2, B is the oscillator frequency, and e = ±1. In this case our r(r) 
rn^B^r^ + + mB (e [2j + 1] — /3), a; = 4, = I /mB, and our leading term reads 



2mB {2k + i + 3/2) + m'^ + mB (e [2j + 1] - /5) 



1/2 



with heigher-order terms identical zeros. Thus, if we take N = 2k + £ ( the harmonic 
oscillator principle quantum number) we come out with the exact Dirac oscillator's closed 
form solution (see equation (35) in [13]) 



E^-m'^ ^ [2mB {N + 3/2) + mB (e [2j + 1] - P)] . 



D. Pure scalar linear potential 

A pure scalar linear potential, i.e., S{r) = Ar and V{r) = 0, is precisely a quark confining 
potential. It has been used by Gunion and Li [5] in Dirac equation to find, numerically, part 
of Dirac J/^' mass spectra. 

Obviously, for this potential equation (20) has to be solved numerically. Then one can 
proceed to obtain the mass spectra for A — 0.137 GeT^^, m — 1.12 GeV, k — — {£ + 1), and 
K = i through the prescription M = 2E. 

In tables 1 and 2 wc report our results for J/^ mass ( in GeV) for k = — (£ + 1) and 
K. = i, respectively. To show the trends of convergence of our results, we report them as 
M(N) — 2E{N), with denoting the number of corrections added to the leading-order 
approximation E^~^\ Our results are also compared with the numerically predicted ones of 
Gunion and Li [5]. Evidently, the accuracy and trend of convergence are satisfactory. 
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E. Funnel-shaped potential 

The funnel-shaped potential is widely used in quarkonium physics. It has both vector 
and components, V{r) and S{r), respectively. 

In Dirac equation, Stepanov and Tutik [4] have used numerical integrations and h- 
expansion formalism without the traditional conversion of Dirac equation into a Schrodinger- 
like form (unlike what we have already done in section II above). They have obtained the 
Charmonium masses for V{r) = —2a/3r and S{r) = br/2, where m = 1.358 GeV , a = 0.39, 
and 6 = 0.21055 GeF^ 

In table 3 we show our results for the Charmonium masses ior k — i and compare them 
with those of numerical integration and /i-expansions of Stepanov and Tutik [4]. They are 
in good agreement and the trend of convergence of our results is also satisfactory. However, 
in table 4 we report the Charmonium masses, for k = — (£+ 1). Therein, we only list our 
results where the mass series and Pade approximants stabihze. 

In KG equation Kobylinsky, Stepanov, and Tutik [3] have used V{r) — —a/r and S{r) — 
br with m = 1.370 GeV , b = 0. 10429 Gel^^ and a = 0.26 to obtain the energies for this 
funnel-shaped potential. They have also used ^-expansions and numerical integrations. In 
table 5 we list our results and compare them with those of /i-expansions, Eji, and numerical 
integrations, Enum, reported in [3]. They are in excellent agreement. 



F. Power-law potential 
An equally mixed scalar and vector power-law potential 

V{r) = S{r) = Ar"" + Bo 

where i/ = 0.1 ( the Martin potential [6] ) and A > is found most successful in describing 
the entire light and heavy meson spectra in the Dirac equation (cf. e.g., Martin 1980, 1981, 
and Jena and Tripati in [6]). Once such potential setting is used in Eq.(6) along with the 
substitution q = r/ g, with 

one gets a simple Schrodinger-type form 



q^ 



n{q) = En{q), 



where 



E = {E-m- 2Bo) {E + m) (2^) 



-2/i/ 



1^/(1^+2) 



(44) 



(45) 



Therefore, one better solve (44) for E (N) and then to find the Dirac quark binding energies 
E from (45). In table 6 we compare PSLET results with those obtained numerically, Enum, 
by Jena and Tripati [6]. The results from the shifted 1/N - expansion technique by Roy and 
Roychoudhury [6], £^i/Ar, are also listed. 
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IV. CONCLUDING REMARKS 



In this work we presented a straightforward extension of our PSLET recipe [7] to solve 
for the eigenvalues of Dirac and KG equations, with Lorentz vector and/or Lorentz scalar 
potentials. We have, again, documented (through tables 1-5) the usefulness of this recipe 
by immediate comparisons with available numerical integration data. 

Nevertheless, one should notice that our results from KG equation are only partially 
better , compared with those from /i-expansions and numerical integrations, than our results 
from Dirac equation. The reason is obviously, and by large, manifested by our approximation 
in equation (8). For Klein-Gordon equation A = in (9) while for Dirac equation A = 1. 

In the process, moreover, there still remain some issues of delicate nature. Namely, 
one can not obtain (using our PSLET above) the exact eigenvalues for V{r) = —Ai/r and 
S{r) = 0, V{r) = and S{r) = —A2/r, or even for V{r) = —Ai/r and S{r) = —A2/r in 
Dirac equation. 

The remedy seems to be feasible in a sort of combination between the current relativistic 
PSLET and a similarity transformation (cf., e.g. ref. [8] and related references therein). 
Preliminary results show that if S{r) — > —A2/r + S'o(r) and V{r) — > —Ai/r + Vo{r) in 
(1) such that So{r) — > and Vo(r) — > as r — > 0, then a similarity transformation 
could accompany our relativistic PSLET to obtain exact results for the generahzed Dirac- 
Goulombic problem and better results for potentials of the type S{r) — —A2/r + So{r) 
and V{r) — —Ai/r + Vo{r). That is, one may carefully follow Mustafa's work [8] to obtain 



+ + U{r) + + - {A2 m(r) + E{r) A^) 



$(r) = £;2$(r). (46) 



where 7 = y k,'^ — Aj + A2 , s = ±1, E{r) = E — K,(r), m(r) = m + So{r), and U{r) — ^ 
as Vo{r) — > and So{r) — > (hence, m(r) — > m and E{r) — > E). Therefore, one 
would replace our in (12), by f = -1/2 7 + s/2 and obtain ( following PSLET 



recipe above) for (48), V{ro) 



m 



2mA2ro, V{r) = -AJr, a; = 2, /3„ = 



l^l - = (^ + 1/2 s / 2 7). In turn, equation (23) yields 



4Q Q-2mA2ro + Al + Am'rl Al- A\ 







to solve for Tq. This would lead to 



m 



Q 



-1/2 



Q 



A; + l/2 + s/2+ 



and 



1 - 



Q 



1/2 



g= (a; + 1/2 + s/2+ y^2^^y 



(/c -I- 1), and 
(47) 

(48) 
(49) 



for = 0, Al 7^ and A\ = 0, A2 ^ 0, respectively. It should be noted that these are the 
well known exact results (cf., e.g.; ref [8]) with the heigher-order terms vanish identically. 

Before we conclude it should be noted that if our results are to be generalized to d- 
dimensions we may incorporate interdimensional degeneracies associated with the isomor- 
phism between the angular momentum £ and dimensionality d (cf., e.g., Mustafa and Odeh 
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(2000) [7]). This would replace our nhy Kd = s{2j + d- 2)/2, where f.,i = f. + {d - 3)/2. In 
this way we reproduce Stepanov and Tutik's [4] and Dong's [9] results in d-dimensions. 

Finally, the above has been a very limited review and a number of other useful and 
novel approaches such as those of Franklin [10, and references therein], Njock et al [11, and 
references therein], • • -etc., have not been touched on. 
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Appendix A: Further algebraic simplifications for PSLET rela- 
tivistic recipe 

Eliminating T-dependence, equation (42) can be simplified into four recursive relations 

to read 



(0)/ 



(50) 



(0), 



and for n > 



(51) 



(52) 



where 



+ E 2Lg (x) Ui:r\x) + x' \U^) (x) + Rt}{x) - v(-\x) 



p=0 



n ^ 



p=Q 



(53) 



(54) 



S^:}{x) = 2kx'-^Gt\{x) + E 2Lg'(x) Gt,7\x) 

p=0 



+x'' 



G^l (x) + Ql:l(x)] + E Lg(x) (Giy) (x) + Qt',7\x)) (55) 



p=0 



Mg(x) = x' pff (x) + E mix) piy^ (x) 

p=0 



(56) 



p=0 



(57) 



Ot}{x) = x*^ [x) + E Lg(x) [x) 



p=0 



(58) 



p=0 



(59) 
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p=0 



p=0 



p=0 



p=0 



fe-i 

p=0 
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Tablel: PSLET results for part of Dirac J/^ spectra (in GeV) for S{r) = Ar, A = 
OmGeV^ V{r) = 0, k = -(^+1), and m = l.UGeV. Where M{N) = 2E{N) with 
denoting the number of corrections added to the leading - order term E^~^^ and 
Mnum are the numerical values reported by Gunion and Li [5] . 



k 


M{N) 


i = 


£= 1 


£ = 2 


i = 3 




Mil) 


3.0919 


3.43078 


3.711960 


3.9581219 




M{2) 


3.0961 


3.43252 


3.712914 


3.9587277 




M(3) 


3.0963 


3.43259 


3.712947 3.9587465 


n 
u 


M(4) 


3.0961 


4.43256 


3.712940 


3.9587436 


M(5) 


3.0961 


3.43256 


3.712939 


3.9587434 




M(14) 3.0961 3.43256 3.712939 3.9587434 






3.103 


3.442 


3.725 


3.973 




Mil) 


4.131 


4.3395 


4.5325 


4.71334 




M(2) 


4.142 


4.3468 


4.5378 


4.71739 




M(3) 


4.148 


4.3502 


4.5401 


4.71905 




M(4) 


4.150 


4.3515 


4.5408 


4.71950 




M(5) 


4.151 


4.3519 


4.5410 


4.71961 


2 


M(6) 


4.152 


4.3520 


4.5411 


4.71965 




M(7) 


4.152 


4.3521 


4.5411 


4.71966 




M(14) 


4.152 


4.3521 


4.5411 


4.71966 




^num 


4.158 


4.36 


4.551 


4.732 
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Table 2: Same as table 1 for k = £. 



k 


M{N) 




£ = 2 


£ = 3 


£ = 4 




M(l) 


3.47090 


3.760125 


4.0111817 


4.2364739 




M(2) 


3.47183 


3.760677 


4.0115488 


4.2367365 




M(3) 


3.47188 


3.760700 


4.0115615 


4.2367444 


(J 


M(4) 


3.47186 


3.760696 


4.0115597 


4.2367435 


M(5) 
: 


3.47186 
: 


3.760695 
: 


4.0115595 
'. 


4.2367435 
: 




M(14) 


3.47186 


3.760695 


4.0115595 


4.2367435 






O /IT 

3.47 


3.757 


4.UU6 


4.23 




M(l) 


3.9570 


4.19083 


4.40310 


4.599080 




M(2) 


3.9624 


4.19451 


4.40578 


4.601126 




M(3) 


3.9640 


4.19537 


4.40631 


4.601482 




M(4) 


3.9644 


4.19557 


4.40642 


4.601542 


1 


M(5) 


3.9646 


4.19561 


4.40644 


4.601552 




M(6) 


3.9646 


4.19563 


4.40644 


4.601554 




M(14) 


3.9646 


4.19563 


4.40644 


4.601554 




^num 


3.965 


4.194 


4.403 


4.597 
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Table 3: PSLET Charmonium masses M{N) = 2E{N) for the funnel-shaped potential, 
V{r) = -2a/3r and S{r) = br/2, with m = 1.358 Gel/, b = 0.21055^6^2^ a = 0.39 
and K = i. The quantum numbers in parentheses are ( k,i). Mnum is the numeri- 
cal integration and M^j is the /i-expansion result ( up to the tenth-order correction) 
reported by Stepanov and Tutik [4]. 



M{N) (0,1) (0,2) (1,1) (1,3) (2,1) (2,3) 



M(l 
M(2 
M(3 
M(4 
M(5 
M(6 
M(7 
M(8 



3.5071 3.8012 3.966 4.3862 4.333 4.6906 

3.5062 3.8007 3.963 4.3857 4.331 4.6908 

3.5056 3.8006 3.961 4.3853 4.329 4.6905 

3.5055 3.8005 3.959 4.3852 4.327 4.6901 

3.5055 3.8005 3.959 4.3851 4.326 4.6899 

3.5055 3.8005 3.958 4.3851 4.325 4.6898 

3.5055 3.8005 3.958 4.3850 4.325 4.6897 

3.5055 3.8005 3.958 4.3850 4.324 4.6897 



M(14) 3.5055 3.8005 3.958 4.3850 4.324 4.6897 

Mn 3.4998 3.7974 3.9501 4.3812 4.316 4.6858 
Mr,um 3.4998 3.7974 3.9499 4.3812 4.315 4.6858 
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Table 4: Same as table 3 for k — — (i+l). Here we report the Charmonium masses where 
the mass series and Pade approximants M^j stabilize. 



k, e M{N) 



M[i,j] k,e M{N) 



M[i,j 



0,0 M(6)=3.0333 M[4,4]=3.0333 1,0 M(7)=3.65 M[5,5]=3.6502 

0,1 M(5)=3.4918 M[2,3]=3.4918 1,1 M(7)=3.946 M[4,4]=3.9462 
0,2 M(4)=3. 7787 M[2,3]=3. 7787 1,2 M(7)=4.1690 M[4,4]=4.1690 
0,3 M(4)=4.0129 M[2,3]=4.0129 2,0 M(9)=4.08 M[6,6]=4.0789 
0,4 M(4)=4.2177 M[2,3]=4.2177 2,1 M(9)=4.314 M[4,5]=4.3139 
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Table 5: KG results for the funnel-shaped potential S{r) = br and V{r) = —a/r, with 
m = 1.370 Gel/, b = 0.10429^61/2, a = 0.26, En represents the results of Kobyhnsky 
et al [3] via /i-expansion (up to the third-order correction), and Enum is the numerical 
integration value reported in [3]. 



E{N) k = OJ. = k = 0,i=l k = 0,i = 2 

E{1) 1.541 1.76167 1.90420 

E{2) 1.535 1.76064 1.90388 

E{3) 1.534 1.76037 1.90380 

E{'i) 1.533 1.76033 1.90379 

E{U) 1.533 1.76033 1.90379 

En 1.536 1.7604 1.9038 

Enum 1.533 1.760 1.904 
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Table 6: PSLET results for E of equation (44) along with the numerically predicted ones 
by Jena and Tripati [6] and the 1/A^-expansion method i?i/Ar by Roy and Roychoudhury 
[6]. in E[N) denotes the number of corrections added to the leading term where 
the series stabilizes. 





E{N) 


E 

-'-'num. 


Ei/N 


0, 


E{2) =1.2358 


1.2364 


1.240 


1,0 


E{7) =1.3347 


1.3347 


1.340 


2,0 


^(4) =1.3922 


1.3923 


1.398 


0, 1 


^(1) =1.3072 


1.3071 


1.309 


1, 1 


^(4) =1.3731 


1.3731 


1.411 


0,2 


£^(1) =1.3540 


1.3544 


1.358 
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